manuscript 



Tunable effective g-factor in InAs nanowire quantum dots 
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We report tunneling spectroscopy measurements of the Zeeman spin splitting in InAs few-electron 
quantum dots. The dots are formed between two InP barriers in InAs nanowires with a wurtzite 
crystal structure which are grown using chemical beam epitaxy. The values of the electron gr-factors 
of the first few electrons entering the dot are found to strongly depend on dot size. They range from 
close to the InAs bulk value in large dots \g*\ = 13 down to \g*\ = 2.3 for the smallest dots. 
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PACS numbers: 73.23.Hk, 73.63.-b, 73.63.Kv, 71.70.Ej 



The spin of an electron in a quantum dot (QD) is one 
of the candidates for a scaleable quantum bit, the funda- 
mental unit in quantum computation and quantum com- 
munication schemes |l| . Experimental realizations are on 
the one hand pursued using top-down approaches. This 
usually involves lateral gate electrodes electrostatically 
confining few or a single electron in a two dimensional 
electron gas close to the surface of a Ga(Al)As based het- 
erostructurc 2J. Such systems offer good tunability and 
controlled coupling of multiple spins has been demon- 
strated [3 . On the other hand, bottom up systems such 
as self assembled QDs Q and carbon nanotubes || are 
expected to scale more easily. Semiconductor nanowires 
have emerged as a promising bottom-up fabricated sys- 
tem for electronic and optical device applications |(| . We 
have recently demonstrated the creation of few-electron 
QDs using InAs nanowire heterostructures with two 
InP barriers. In the following we set out to investigate 
the spin properties of the first few orbital levels of these 
QDs. 

We utilize transport spectroscopy to measure the Zee- 
man splitting of the energy levels as a function of mag- 
netic field and thereby determine the effective electron 
g-factor (g*). The g-factor of bulk InAs, which crystal- 
izes in the zinc-blende (ZB) structure, has been found 
to be g* = —14.7 |8|]. However, InAs nanowires can ex- 
hibit both zinc-blende and wurtzite (WZ) type crystal 
structure depending on diameter and growth condi- 
tions and so far very little is known about band parame- 
ters in WZ InAs. In low-dimensional semiconductor het- 
erostructures the g factor depends critically on system 
size and dimensionality [Tcj| . We show that varying the 
size of our nanowire dots allows us to tune g* from a value 
close to the InAs bulk value down to \g*\ = 2.3 ±0.3. The 
possibility to have multiple dots along a nanowire, each 
with a different g-factor, makes such systems interest- 
ing candidates for realizations of individually addressable 
spin qubits. 

Using chemical beam epitaxy InAs nanowires contain- 
ing QDs were grown catalytically from Au nanoparticles 
deposited on a <111>B InAs substrate [HI E|- The 
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FIG. 1: (a) STEM image of a nanowire together with a 
schematic of the InAs QD of diameter D. (b) High resolu- 
tion STEM image of a QD from the same growth as S#4 and 
S#6. The dot length L =12 nm is deduced from the length 
of the InAs segment (bright) between two InP barriers (dark) 
having widths wi and W2- (c) Coulomb blockade oscillations 
in Lsd as a function of V g at B =0 T and a bias Vsd = 50 /xV 
for S#4. Each peak corresponds to adding a single electron to 
the dot. (d) Dot conductance as a function of magnetic field 
and gate voltage. The peaks shift as a function of B because 
of the Zeeman effect and small B induced orbital shifts, (e) 
Peak separation AVg as a function of B showing the Zeeman 
splitting for the first electrons added to the dot. 
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nanowires typically grow perpendicular to the substrate 
and high resolution scanning transmission electron mi- 
croscope (STEM) images indicate that most of them 
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TABLE I: Summary of important parameters for the different 
devices measured, including the measured effective g-factors 
averaged over data extracted from peak position and excited 
state spectroscopy. 



show a nearly perfect WZ lattice with very few twinning 
boundaries or stacking faults. The QDs were defined be- 
tween two InP tunnel barriers within the InAs nanowires. 
This yields a confining potential for electrons set by the 
600 mcV conduction band offset between InAs and InP in 
the growth direction ^1 ■ in the lateral direction the side 
facets of the wire form a hexagonal cross section with pre- 
sumably hard wall conditions [see Fig. [Ja)]. The wires 
studied here had diameters D in the range of 50 — 70 nm. 
Five growths were made with wires containing QDs of 
different lengths. The dot length L and the thicknesses 
w\ and W2 of the InP tunnel barriers were measured for 
each growth using STEM images such as the one shown 
in Fig. n^b). A small variation in D within the same 
growth (given by the size distribution of the gold cata- 
lyst particles) leads to slightly different growth rates [TT| . 
The STEM images therefore allow us to establish a rela- 
tion between D and L,w\,W2 for each growth. We then 
measure the diameter D of each nanowire in the actual 
device with a scanning electron microscope, which allows 
us to give estimates for all the size parameters accurate 
to within about 15% (see Table P). 

For transport measurements the samples were pre- 
pared by depositing wires on a degenerately doped Si sub- 
strate with a 100 nm SiC>2 layer acting as the gate dielec- 
tric for the Si back gate. After locating individual wires, 
Ni/Au source and drain contacts were defined using elec- 
tron beam lithography and metal evaporation. Device 
S#8 was fabricated as a reference from a sixth growth 
of homogenous InAs nanowires without QDs. Here, local 
gates were used to define large QDs {L « 200 nm) [l3| • 

Electrical measurements were performed under sym- 
metric bias conditions in a He-3 cryostat with a base 
temperature of 250 mK. The heterostructure dot lengths 
range from L = 8 — 20 nm which is much smaller than 
D. This leads to strong quantum confinement along the 
wire, so that only the lowest quantum state is occupied 
in the growth direction. A magnetic field B is applied 
perpendicular to the wire. 

Figure ^c) shows Coulomb blockade (CB) oscilla- 
tions in the current I$d through the QD in S#4 with 
L =12 nm. At zero applied gate voltage the QD is empty 
(N = 0). For V g — 686 mV the lowest orbital level of the 



QD aligns with the chemical potential in the source and 
drain contacts giving rise to a peak in I$d ■ Upon adding 
more electrons the Coulomb peak positions on the V g axis 
directly reveal the addition spectrum of the QD for which 
a clear shell structure is observed Q ■ Up to at least the 
10th electron the shell structure is solely determined by 
the lateral confinement. For higher electron numbers the 
second subband gets occupied for the largest heterostruc- 
ture dot with L — 20 nm. The enhanced peak spacing 
between the second and third peak corresponds to a filled 
first shell i.e. the lowest orbital containing both a spin up 
and a spin down electron. The spin degeneracy is lifted 
at finite B due to the Zeeman splitting AE Z = g*nsB 
where /is is the Bohr magneton. The Zeeman shift of 
the energy levels manifests itself directly as a change 
in the CB-peak positions 01 • Figure ^d) shows the 
CB peaks as a function of B where peaks reflecting the 
charging of the same orbital level with spin up and spin 
down electrons move apart. Such spin pairing, indicat- 
ing alternating spin filling for the first few orbital levels, 
is observed for all our dots except the two largest ones 
(S#7 and S#8) in strong magnetic fields. In addition to 
the Zeeman effect the peak positions are also shifted due 
to a small diamagnetic shift. Such contributions can be 
eliminated by taking the difference between two peak po- 
sitions for which the same orbital is filled i.e. spin paired 
peaks neighbouring a gap with odd electron number. The 
resulting peak spacing AV g is proportional to the addi- 
tion energy U(N) which can be separated into three parts 
U(N) = E mt + AE(N) + AE Z . Here E int is the interac- 
tion contribution, which within the constant interaction 
modelQ is simply the charging energy E int — e 2 /Cdot, 
AE(N) is the single particle level spacing and AE Z is 
the Zeeman splitting. The capacitive lever arm between 
V g and dot energy can be extracted for each CB peak 
from measurements of Coulomb diamonds and is found 
to be in the range of a = Cdot/C g = 0.05 — 0.25 de- 
pending mainly on dot size. The effective g-i actor is now 
given from the slope of the curves obtained by calculat- 
ing AVg(B) as shown in Fig.QJe) for three different dots. 
The topmost graph gives the first five peak spacings for 
dot S#4, where the spin pairs are given by the positive 
slopes corresponding to odd numbered gaps [see legend 
in Fig^e)]. Fits to the measured data for these spin 
pairs yield \g*\ — 4.4 ± 0.5. The same type of analysis 
was done for all eight samples giving g-factors between 
= 2.3 ± 0.3 for a dot with L = 8 nm (S#l) and 
\g*\ = 8.9 ± 1.0 for a dot with L = 20 nm (S#7). For 
the latter case the large spin splitting leads to Zeeman 
induced crossings, which explains the deviations from the 
linear behaviour of AV g . In addition, we were not able to 
measure the first two diamonds with sufficient accuracy 
in this sample, which is why the corresponding lines are 
missing in the data. 

We also determine the Zeeman splitting directly from 
the excited state spectrum obtained in Coulomb diamond 
measurements. Fig. |2Ja) shows the differential conduc- 
tance G as a function of Vg and Vsd in grayscale at 



3 



(a) 



0.9 



N S B=10T 


1 1 










* i *7 — i — 1 1 


N=0 |SL- 

1 u 




j : 
> ! 




i 




N=0 


-0.01 c 


) 0.01 



> 
> 

0.8 
0.7 
0.6 

(b)i 2 

p 8 

S 4 


-u.u I 

FIG. 2: (a) Differential conductance as a function of gate 
voltage and applied bias for S#3 at B = 10 T. The Zeeman 
split excited states are seen as lines running parallel to the 
diamond edges as indicated by the arrows, (b) Differential 
conductance at a gate voltage of 620 mV as a function of 
bias and magnetic field where the linear Zeeman splitting is 
directly observed. 



B = 10 T again for S#4. For the lowest open diamond 
the dot is empty (N = 0). The borderlines of this dia- 
mond mark the situation where either the source or drain 
chemical potential is aligned to the first dot level be- 
ing charged with a single spin up electron. The Zeeman 
splitting is then observed as an excited state line mov- 
ing parallel to the borderlines but shifted to larger gate 
volt ages |2J (see arrows). The same is true for excitations 
below and parallel to the borderlines of the N = 1 to 
N = 2 transition. The corresponding transitions are not 
observed within the N — 2 diamond since the first shell 
can accommodate two spins only and the total spin of 
the dot is expected to be zero in that situation. The 
lines outside the diamonds that have a different slope 
than the diamond edges and show negative differential 
conductance (white lines) are due to a modulation of 
transport through the dot by states in the quasi one- 
dimensional source and drain contacts. In order to make 
the difference between these lead states and excited state 
lines more clear, Fig.|2b) shows G as a function of B and 
source-drain bias for fixed Vg — 620 mV [dashed white 
horizontal line in Fig.|^a)] . Here the linear Zeeman split- 
ting is directly observed as a function of B. The lead 
states appear as a fast modulation of the conductance 
and to a lesser degree also influence the position of the 
diamond borderlines. The splitting at B = 10 T is out- 
lined by vertical dashed lines. Measuring the splitting as 
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FIG. 3: Dependence of dot length on the effective g-factor. 
The top curve is the theoretical values obtained from Eq. Q 
and hard wall confinement. 



a function of B in Fig. Hfb) we deduce \g* | = 4.3 ±0.5 in 
good agreement with \g*\ = 4.4 ± 0.5 found above from 
peak separations. 

The effective g-factors for 8 dots of different lengths 
were deduced in the same way as above and are plotted 
as a function of L in Fig. EI Here, L denotes the extent of 
the dot in the strongest confinement direction. For S#8, 
which was the dot defined using local gate electrodes, we 
therefore replaced L by the diameter D of the nanowire 
since the separation of the gate electrodes was more than 
200 nm. We find that in this case g* approaches nearly 
the bulk value found for ZB InAs. For the heterostructure 
dots with small L, g* is strongly suppressed. 

Generally, g-factors in confined semiconductor systems 
have been found to differ si gnifi cantly from the corre- 
sponding bulk values 0, ITH " ll6t Il7j . In measurements 
on In As-based two-dimensional electron systems, values 
from g* — — 1 to g* = — 13 have been found experimen- 
tally depending on magnetic field direction |l 8l. a uan- 
tum well width [15j and barrier composition llaJlTI and 
gate tunability of g* has been demonstrated [171 Il9|. In- 
vestigations on self-assembled InAs QDs using magneto- 
tunneling and capacitance spectroscopy have revealed 
anisotropic g-factors with values from +0.5 to +1.6 [2(], 
21]. For WZ crystals both effective mass and g-factor 
are expected to be anisotropic and theoretical expression 
can be derived within k ■ p-theory[23, Hj| . However, to 
the best of our knowledge, band parameters for WZ InAs 
are not known. We therefore, start the discussion of our 
results with a comparison to a simple expression derived 
for cubic III / V-compounds . 
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Here g* depends on the bandgap E g of the semicon- 
ductor material and on the spin-orbit splitting A so in 
the valence band. P 2 is a band parameter related to 
the interband momentum matrix elements. The sup- 
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pression of g* can be understood qualitatively in a sim- 
ple picture in which the dot confinement leads to an en- 
hanced gap between electron and hole states. We there- 
fore assume that E g = E h g wlk + E\ where E\ is the con- 
finement energy for the lowest subband of the dot cal- 
culated using a hard wall confinement potential with 
width L. For WZ semiconductors the gap is typ ically 
slightly larger than for the same ZB material 24] . Ac- 
cordingly, we use 

E buik = o 450 eV a value which is 10% 
larger than the ZB value. Beyond that, we rely on the 
band parameters for ZB InAs which we expect to be a 
good first approximation[2jJ 0] ■ In Fig. the solid line 
shows a curve derived from Eq. (JTJ with P 2 = 21.5 eV, 
A so = 0.39 eV and a constant effective electron mass of 
0.026[2^. While this model reproduces the general trend 
of the measured data points and fits reasonably well for 
the larger dots, there is a much stronger suppression of 
g* with decreasing L in the experiment. 

We have considered the reduction in g* for small dot 
sizes to be due to penetration of the wavefunction into 
the InP barriers which have a bulk g- factor of +1.40. 
However, in the experiment we observe no measurable 
change in g* with electron number. For large B the over- 
all conductance is observed to strongly decrease due to 
a stronger magnetic confinement which results in weaker 
coupling but no observable change in g*. In Table [i] we 
also give values for the Coulomb peak conductance G at 
B = T and small bias (average over first five peaks) 
which is clearly correlated to the barrier thickness. This 
demonstrates tunability of the coupling while keeping L 
fixed. We do not find, however, that this is correlated 
with a change in g* (see e.g. S#3 and S#4). This leads 
us to conclude that wavefunction overlap with the barrier 



material is unlikely to be the cause of the more strongly 
reduced values of g* for small L in the experiment. A 
more accurate model is necessary, taking into account 
both the confinement and the WZ crystal structure in 
order to explain the observed data. 

In conclusion we have measured the effective electron 
g-factor of InAs QDs defined by InP double barriers in 
semiconductor nanowires with a WZ crystal structure as 
a function of dot size. We observe a strong suppression 
of g* for small dot length L down to \g*\ — 2.3 for S#l 
(L = 8 nm). While this trend can be understood qual- 
itatively due to an increase of the confinement energy 
a good theoretical model for WZ InAs is missing. We 
hope that our experiments will spur the development of 
more quantitative models for these type of InAs QDs. 
It was further shown that changing the dot length L al- 
lows us to design QDs along a nanowire with different 
specific spin splittings in a constant magnetic field. Fu- 
ture development of graded heterostructures is ex pect ed 
to allow individual gate tunable spin splittings p^.ll9| in 
a series of dots along a nanowire, which makes nanowire 
QDs containing a single electron spin interesting systems 
for the realization of qubits. 

Acknowledgments 

This work was supported by the Swedish Founda- 
tion for Strategic Research (SSF), the Swedish Research 
Council (VR), the Knut and Alice Wallenberg Founda- 
tion (KAW), the Office of Naval Research (ONR), and 
the Swiss Science Foundation (Schweizerischer Nation- 
alfonds). Helpful discussions with C. Thelander and S. 
Ulloa are gratefully acknowledged. 



[1] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

[2] R. Hanson, B. Witkamp, L. M. K. Vandersypen, L. H. W. 
van Beveren, J. M. Elzerman, and L. P. Kouwenhoven, 
Phys. Rev. Lett. 91, 196802 (2003). 

[3] N. J. Craig, J. M. Taylor, E. A. Lester, C. M. Mar- 
cus, M. P. Hanson, and A. C. Gossard, Science 304, 565 
(2004). 

[4] M. Kroutvar, Y. Ducommun, D. Heiss, M. Bichler, 
D. Schuh, G. Abstreiter, and J. J. Finley, Nature 432, 
81 (2004). 

[5] M. Biercuk, S. Garaj, N. Mason, J. Chow, and C. Marcus, 

cond-mat p. 0502634 (2005). 
[6] L. Samuelson, Materials Today p. 22 (2003). 
[7] M. T. Bjork, C. Thelander, A. E. Hansen, L. E. Jensen, 

M. W. Larsson, L. R. Wallenberg, and L. Samuelson, 

Nano Lett. 4, 1621 (2004). 
[8] Landolt-Bornstein, Physics of group IV elements and III- 

V compounds, vol. 17 (Springer- Verlag, 1982). 
[9] M. Koguchi, H. Kakibayashi, M. Yazawa, K. Hiruma, 

and T. Katsuyama, Jpn. J. Appl. Phys. 31, 2061 (1992). 
[10] A. A. Kiselev, E. L. Ivchenko, and U. Rossler, Phys. Rev. 

B 58, 16353 (1998). 



[11] L. E. Jensen, M. T. Bjork, S. Jeppesen, A. I. Persson, 

B. J. Ohlsson, and L. Samuelson, Nano Lett. 4, 1961 
(2004). 

[12] M. T. Bjork, B. J. Ohlsson, T. Sass, A. I. Persson, 

C. Thelander, M. H. Magnusson, K. Deppert, L. R. Wal- 
lenberg, and L. Samuelson, Appl. Phys. Lett. 80, 1058 
(2002). 

[13] C. Fasth, A. Fuhrer, M. T. Bjork, and L. Samuelson, 

Nano Letters (2005). 
[14] L. P. Kouwenhoven, D. G. Austing, and S. Tarucha, Rep. 

Prog. Phys. 64, 701 (2001). 
[15] B. Kowalski, P. Omling, B. K. Meyer, D. M. Hofmann, 

C. Wetzel, V. Harle, F. Scholz, and P. Sobkowicz, Phys. 

Rev. B 49, R14786 (1994). 
[16] M. Dobers, K. v. Klitzing, and G. Weimann, Phys. Rev. 

B 38, 5453 (1988). 
[17] J. Nitta, Y. Lin, T. Akazaki, and T. Koga, Appl. Phys. 

Lett. 83, 4565 (2003). 
[18] S. Brosig, K. Ensslin, A. G. Jansen, C. Nguyen, B. Brar, 

M. Thomas, and H. Kroemer, Phys. Rev. B 61, 13045 

(2000). 

[19] G. Salis, Y. Kato, K. Ensslin, D. C. Driscoll, A. C. Gos- 
sard, and D. D. Awschalom, Nature 414, 619 (2001). 



■5 



[20] A. S. G. Thornton, T. Ihn, P. C. Main, L. Eaves, and 
M. Henini, Appl. Phys. Lett. 73, 354 (1998). 

[21] G. Medeiros-Ribeiro, M. V. B. Pinheiro, V. L. Pimentel, 
and E. Marega, Appl. Phys. Lett. 80, 4229 (2002). 

[22] C. Hermann and C. Weisbuch, Phys. Rev. B 15, 823 
(1977). 

[23] M. W. Bayerl, M. S. Brandt, T. Graf, O. Ambacher, 



J. A. Majewski, M. Stutzmann, D. J. As, and K. Lischka, 
Phys. Rev. B 63, 165204 (2001). 
[24] M. Murayama and T. Nakayama, Phys. Rev. B 49, 4710 
(1994). 

[25] I. Vurgaftman, J. R. Meyer, and L. R. Ram-Mohan, J. 
Appl. Phys. 89, 5815 (2001). 



